Spin stiffness and quantum fluctuations in C-type and A-type antiferromagnets 
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We present a systematic study of quantum fluctuations in the C-type and A-type antiferromag- 
netic (AF) phases in cubic lattices and in bilayer systems. Using the linear spin-wave theory, we show 
that the spin stiffness and the quantum corrections to the order parameter and energy obtained for 
C-AF and A-AF phases decrease with the increasing number of ferromagnetic bonds. Therefore, the 
quantum spin effects in LaMnC>3 and in LaVC>3 are rather small, suggesting the magnetic moments 
of ~ 3.91 and ~ 1.89/is, respectively. They cannot explain the strong reduction of the magnetic 
order parameter observed in cubic vanadates. [Published in Phys. Rev. B 66, 094431 (2002)] 
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The undoped transition metal oxides are characterized 
by large local Coulomb interactions cx U, which lead to 
several fascinating phenomena such as high-temperature 
superconductivity and "colossal magnetoresistance" - 1 
When the Coulomb interactions dominate over the ki- 
netic energy, the charge fluctuations are quenched and 
the magnetic properties follow from the effective low- 
energy superexchange interactions. In some of these 
systems the orbital degrees of freedom play a role due 
to the partial filling of (almost) degenerate d orbitals, 
the superexchange interactions are strongly frustrated, 2 ' 3 
and the quantum effects are enhanced. 4 These interac- 
tions together with the Jahn- Teller effect may induce 
the orbital ordering below a structural transition and 
break the cubic symmetry of the perovskite lattice. In 
such systems, although the crystalographic directions in 
a three-dimensional (3D) cubic lattice are a priori equiv- 
alent, one finds magnetic interactions not only of differ- 
ent value, but even of different sign, stabilizing C-type 
or A-type antiferromagnetic (AF) phases. 3 

One of the best known examples of the non-cubic mag- 
netic interactions in a perovskite system is the A-type AF 
order observed in LaMn03, 5 with ferromagnetic (FM) 
superexchange within (a, b) planes (J a b) and AF interac- 
tions along the c axis ( J c ), or in KCUF3, an almost perfect 
one-dimensional (ID) Hciscnberg antifcrromagnet. 6 In 
both above cases the magnetic ordering is supported by 
the orbital ordering which is induced either by the Jahn- 
Teller effect, 7 ' 8 or by the superexchange interactions. 9 ' 10 
Recently it was suggested that the latter contribution 
dominates, 10 but this issue is still controversial and has to 
be clarified by future studies. The spin waves in LaMnOa 
have been investigated in great detail and it was estab- 
lished that the AF interactions J c are weaker than the 
FM J a b ones, 9 in good agreement with the experimental 
data. 11 

An inverse situation with respect to manganites and 
cuprates, with AF interactions within (a, b) planes co- 
existing with FM superexchange along the c-axis, is en- 
countered in the so-called C-AF phase, observed in cu- 
bic vanadates: in LaV03 below the Neel temperature 
Tjv — 143 K, 12 and in YVO3 at intermediate temper- 



atures 77 < T < 116 K. 13 Finally, G-type AF or- 
der, with cubic symmetry and the magnetic order pa- 
rameter staggered in all three directions, is found in 
CaMnC>3, LaTiC>3, and also in the low-temperature phase 
of YVO3. 13 Particularly this last example shows that the 
type of magnetic order observed in transition metal ox- 
ides may be triggered by a delicate balance of magnetic 
interactions induced by the orbital ordering. 

Other examples of non-cubic antiferromagnets are 
found in bilayer systems, where the effective dimension- 
ality is reduced and the cubic symmetry is explicitly 
broken by geometry, even when the exchange interac- 
tions are AF and identical on the bonds in different 
crystallographic directions. The G-AF structure is re- 
alized in YBa2Cu306+x, 14 while the bilayer manganites 
La2-2xSri + 2^Mn207 show interesting dependence of the 
magnetic order on the doping, with FM structure for 
x ~ 0.4 (Ref. 15) and A-AF structure for x ~ 0.48, 16 
separated by a C-AF phase at intermediate doping. 17 

The renewed interest in the magnetic properties of 
transition metal oxides motivates a systematic study of 
quantum fluctuations in different AF structures. The 
spin-wave theory was introduced long ago, 18 and high 
accuracy results of the 1/5 expansion were presented for 
a two-dimensional (2D) Heisenberg antiferromagnet. 19,20 
More recently the properties of 2D dimerized models were 
studied. 21 However, we are not aware of any systematic 
investigation of the spin stiffness and the quantum ef- 
fects in non-cubic systems with different strengths and 
signs of exchange interactions. In this paper we con- 
sider the quantum fluctuations in the AF phases real- 
ized in transition metal perovskites, and compare them 
with those known for the 2D square and 3D cubic lat- 
tice. Further motivation comes from the experimental 
data for numerous systems which could be properly un- 
derstood only when spin effects would be extracted from 
the full quantum problem. Therefore, we neglect the or- 
bital fluctuations, and study the spin quantum effects 
alone in the structures stabilized by the superexchange 
fixed by a rigid orbital background. 

We consider the (effective) 3D Heisenberg model with 
nearest-neighbor superexchange interactions, J a b in {a, b) 



1 



planes and J c along c axis, between spins S in G-AF, G- 
AF and A-AF phases, given by 

'Hs.D = Jab ^ Sj • Sj + J c ^ S, • Sj. (1) 

W>ll(<»>6) W)l|c 

In cubic crystals the G-AF phase is obtained for J ab > 
and J c > 0, while FM interactions stabilize either G-AF 
phase if J c < 0, or the A-AF phase if J ab < 0. Simi- 
lar phases are possible in the bilayer structures, with the 
second sum involving only the interlayer bonds. We will 
investigate the range of parameters with the interlayer 
coupling | J c | < 2J ao , where the ground state is ordered. 
In contrast, for J c / J ab > 2.55 the magnetic properties are 
dominated by the singlets which form on interlayer bonds 
(ij) || c axis, the long-range order is lost, and the spin- 
wave expansion does not apply. 22 As a reference system, 
we also performed calculations for the 2D square lattice 
with different interactions along the bonds parallel to a 
and b axis, respectively, described by 
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We assume J a > and study both the G-AF (J b > 0) 
and G-AF ( J b < 0) phase. 

We determined the dispersion of spin waves in the spin 
models given by Eqs. (1) and (2), the quantum correc- 
tions to the magnetic order parameter (S z ), and to the 
ground state energy E = (H^d) (E = (H2d)) using the 
linear spin-wave theory (LSWT) in the leading 1/5 or- 
der. First we transform the spin operators at the sites 
j e B sublattice: Sf -» Sf, S] -» -S], which removes 
spins down (»S|)o = — S of the Neel state. Next we intro- 
duce bosonic operators {a\, c^}, and use the lowest order 
of the Holstcin-Primakoff transformation: 23 



St 



2Sa, 



2Sa\. (3) 



The excitations are derived from the equations of motion 
for the energy-dependent Green's functions, 24,25 

E(^i\a])) = ^S ij + (([a i ,n 3 D}\a]}}. (4) 

After the Fourier transformation the spin waves are found 
by a Bogoliubov transformation which diagonalizes a 2 x 2 
dynamical matrix at each k point. The positive magnon 
energies for the 3D G-AF, G-AF and A-AF phase are: 

w G (k) = 25{(2J ab + J c f 

-[2J abl+ (k) + J^k)] 2 } 1 / 2 , (5) 
w c (k) = 25{[2J ah + | J c |(l - 7z (k))] 2 

-^(k)} 1 / 2 , (6) 
w A (k) = 2S{[2\ J ab \{\ - 7+ (k)) + J c ] 2 



where 7+(k) = i (cos k x + cos k y ) and 7z(k) = cosfc z 
are k-dependcnt structure factors. The corresponding 
magnon energies for the bilayer systems are: 

4(k) - S{(4J ab + J c ) 2 - [4J ah7+ (k) + AJJ 2 } 1 / 2 , (8) 
4(k) = S{[4J ab + | J c |(l - A)] 2 - lej^^k)} 1 / 2 , (9) 
4(k) = 5{[4| J o6 |(l - 7+ (k)) + J c ] 2 - A 2 J 2 } 1 / 2 , (10) 

where A = ±1 corresponds to k z — 0, n, i.e., to the sym- 
metric and antisymmetric interlayer modes, respectively. 
The magnon energies for the G-AF and G-AF phase in 
the 2D case are: 

c4 D (k) = 2S{(J a + Jb) 2 + [Ja~fx(k) + J b ~/y(k)} 2 } 1 / 2 , (11) 

4 c (k)^25{[J Q + |J b |(l- 7 ,(k))] 2 -J 2 7 2 (k)} 1 / 2 , (12) 

with 7;c (k) = cosfcz, and 7y (k) = cosk y . 

The size of quantum fluctuation corrections to the clas- 
sical order parameter (S'f)o = S, and the intersite spin 
fluctuations oc (Sf Sf) which modify the energy, deter- 
mine the stability of the classical phases. One finds that 
the order parameter for i e A sublattice is reduced by 
local quantum fluctuations, 



{S?) = S-(SrS+) = S-6(S*), 



(13) 



and the local correlation function (S i Sf) is determined 



by the Green's function ((» k |» k )) by, 

i r+oo 
i. J— oo 
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exp(f3ui) — 1 ' 
where j3 = l/fc^T, N is the number of sites, and 
J(k,w) = 2 Im((a k |aj c » u; _ i£ 



(14) 



(15) 



is the spectral density of spin excitations, with the op- 
erator a k being a Fourier tranform of a i . The quan- 
tum corrections are found by taking T — > limit of Eq. 
(14), where the averages (Sf Sf) are determined from the 
spectral weights *4^(k) in Eq. (15) of the excitations at 
negative frequencies, and therefore, 



(16) 



k u<0 



Thus, using Eqs. (15) and (16) one obtains the quantum 
corrections S(S Z ) for the 3D phases: 



HS')c = sJ- 

6(S')c = S 



d 3 k 2J ab + J C _1 
2' 



) 3 w G (k) 

d 3 k 2J afc + |J e |[l~ 7z (k)] 1 

(2tt) 3 Lu c (k) 2' 
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(17) 
(18) 
(19) 
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with magnon dispersions given by Eqs. (5)— (7). The 
quantum corrections for the bilayer structures and for 
the square lattice can be obtained in a similar way - 
they involve 2D integrations over the respective excita- 
tion spectra, and, in addition, the summation over the 
symmetric and antisymmetric modes for the bilayer sys- 
tem. 

The quantum corrections 8E to the ground state en- 
ergy per site, E = Eq — 5E, where Eq is the energy found 
using the classical state, were obtained by integrating the 
excitation spectra wmM, with M = A,C,G. For in- 
stance, for 3D systems one finds E a = — (2| J a b\ + \ J C \)S 2 , 
and 



6E = (2\J ab \ + | J, 



d 3 k 



wjw(k). 



(20) 



In order to compare the systems of different dimensional- 
ity we use below the relative energy correction S(5E/Eq) 
as a measure of the quantum correction to the ground 
state energy. For the 3D AF systems with \J a b\ — \J C \ — 
J (or for the 2D systems with \J a \ — \Jt\ — J) the rela- 
tive energy correction SE/E can be expressed by, 



SE 
Eq 



J_ 

zS 1 



(21) 



where z is the number of nearest neighbors, and the co- 
efficient 7 > provides a measure of quantum effects. It 
depends on the system dimensionality, on lattice type, 
and on the type of AF order. Therefore, we write the 
ground state energy as follows: 23 



E = --zJS 2 
2 



(l • - 



zS) 



(22) 



Consider first a cubic crystal with equal strength of AF 
and FM superexchange interactions (|J &| = \Jc\ = J)- 
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FIG. 1. The spin-wave dispersions as obtained within the 
LSWT in the G-AF, C-AF and A-AF phases along the main 
directions of the 3D BZ for a simple cubic lattice, with 
equal AF and FM exchange interactions (|J a f>| = \J C \ = J)- 
We have used the standard labels for high-symmetry points: 
r = (0,0,0), X = (tt,0,0), M = (7T,7T,0), R = (7T,7r,7r), and 
Z = (0,0,7r). 



The spin-wave excitations with positive energies are pre- 
sented in Fig. 1 along the high-symmetry directions 
of the 3D Brillouin zone (BZ). One finds an expected 
Goldstonc mode with uomQ*) = at R = (n , ir , n) , 
M = (7r,7r,0), and Z = (0,0, n) point for G-AF, C- 
AF and A-AY phase, respectively. It corresponds to the 
vector in the reciprocal space which couples two spin ex- 
citations in a given AF structure. The spin excitations 
are linear near the Goldstone mode, u>mQ<-) — Dk, and 
close to the V point in each AF structure, and the spin- 
wave stiffness constant D is isotropic, being the same for 
any cubic direction, in spite of the different signs of su- 
perexchange constants. Note also that the width of the 
magnon dispersion increases with the increasing number 
of FM bonds, and would be the largest and equal to 24 J 
for the cubic FM phase, with equal FM exchange interac- 
tions in all three directions (not shown) . Increasing over- 
all dispersion UM(k) in Eq. (20) reduces the quantum 
correction 5E when the number of FM bonds increases. 

We compare the quantum effects given by Eqs. (16) 
and (21) for different systems listed in Table I. The val- 
ues obtained for 2D and 3D G-AF phases reproduce the 
known results given by Mattis. 23 Both 5(S Z ) and SE are 
reduced when some of the bonds are FM. When the di- 
mensionality is fixed, the quantum corrections decrease 
in a systematic way with the increasing number of FM 
interactions, from the G-AF to A-AF order. 26 A some- 
what counterintuitive result is that at the same time the 
stiffness constant D is enhanced only by the AF inter- 
actions, while it decreases with the increasing number of 
FM bonds for all systems: the 3D cubic (see Fig. 1), 
2D square, as well as for the bilayer structure (Table 
I). Thus, the magnon energy increases slower with k 
near the T point by a factor l/y/3 in A-AF than in G- 
AF phases. Further, for a given magnetic structure the 
quantum corrections S(S Z ) and S(5E/E ) gradually de- 
crease when the dimensionality of the system increases 
from a 2D square lattice thoughout a bilayer system to a 
3D cubic lattice (Table I). In the bilayer case an isotropic 

TABLE I. The quantum corrections to the magnetic order 
parameter 5{S Z ) and to the ground state energy SE, and the 
stiffness constant D, as obtained in the LSWT for the phases 
with isotropic low-energy spin excitations ojM(k) ~ Dk, with 
M = G,C,A. Parameters: \J ab \ = \J C \ = J for the 2D and 
3D AF phases, and \J a b\ = | J c |/2 = J in the bilayer systems. 
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FIG. 2. Quantum corrections to: (a) the magnetic order 
parameter 8{S Z ), and (b) energy S(6E/Eo), obtained in the 
LSWT as functions of a = \J b \/J a for the G-AF and C-AF 
2D phase. Part (c) shows the stiffness constants: Db and D a 
for G-AF phase (solid and dot-dashed line), and D b for G-AF 
(long-dashed line), normalized to the value Dq D — 2\f2JS 
obtained for the G-AF phase with J a = J b . The LSWT re- 
sult for S(SE/Eo) in the ID model is shown in (b) by dotted 
line. 



FIG. 3. Quantum corrections as in Fig. 2, but for the 3D 
AF phases, and with a = \J C \/J a b {a = \J ab \/Jc) for the G- 
and G-AF phase (A-AF phase). Dotted lines show the values 
of 8(S Z ) (a) and SSE/E Q (b), obtained within the LSWT for 
the isotropic 2D square lattice. Part (c) shows the stiffness 
constants D 1 normalized to Dq D = 2 V / 3JS' obtained for the 
G-AF phase with cubic symmetry of exchange interactions: 
D c and D ab (solid and dot-dashed line) for G-AF phase; D c 
for G-AF phase (long-dashed line), and D ab for A-AF phase 
(dashed line). 



stiffness constant is obtained only when the interlayer 
superexchange interaction is larger by a factor of two 
than the respective intralayer interactions, i.e., |J a &| = 
\J c \/2 — J, which simulates the missing second neighbor 
of each atom along the c axis. 

The long-range order, with (S z ) > (S(S Z ) < 0.5 for 
5=1/2) [see Fig. 2(a)], is stabilized within the LSWT in 
quasi-lD systems realized on a 2D square lattice already 
by relatively small interchain couplings, \ Jb\/J a — 0.15. 
Qualitatively the values of 5(S Z ) and SE/E [Fig. 2(b)] 
behave in a similar way; they both decrease for the C-AF 
phase, while the quantum corrections reach a minimum 
for an isotropic 2D antiferromagnet at J b = J a in the 



G-AF phase, and next increase again for a > 1 when a 
stronger AF interaction occurs along the b axis. 

The magnon stiffness increases when the (either AF or 
FM) interactions along the b axis |Jb| increase, and one 
finds D b /Dl D ~ for small a = \ J b \/J a [Fig. 2(c)], 

with Dq D — 2V2J5' obtained in a G-AF phase with cu- 
bic symmetry. However, the behavior of the G-AF and 
C-AF phase is qualitatively different. First of all, the 
value of D a increases with J b in the G-AF phase, while it 
remains constant {D a /Dl D = 1/V%) in the C-AF phase. 
Second, except for the asymptotic regime of a < 0.3, the 
values of D b are considerably lower in the C-AF phase 
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FIG. 4. Quantum corrections (a) to the magnetic order 
parameter 8{S Z ) (solid line), and to the ground state energy 
SSE/Eo (dashed line), and (b) the stiffness constants D c and 
D ab (solid and dashed line) for an AF bilayer system with 
J a b > 0, as functions of the interlayer coupling J c /J a b, in 
units of D|j = 2^3 J S (see Table I). 

than those in the G-AF phase. 

Next we consider C-AF and A-AF order in 3D systems 
with increasing FM interactions and compare them with 
the G-AF phase (Fig. 3). Both 6(S Z ) and 6E/E de- 
crease rapidly in G-AF phase with increasing a — J c / J ab 
in the regime of < a < 0.5, while for a > 0.5 the 
quantum effects are reduced stronger by the FM bonds 
which make the magnetic order in the G-AF phase more 
classical than in the G-AF phase. The values of S(S Z ) 
and 5E/E increase again for a > 1 in the G-AF phase 
as the system starts to approach a quasi-lD regime with 
enhanced quantum fluctuations. In the A-AF phase the 
strongest quantum fluctuations are obtained for small 
ol = | J a b | /Jc ) where the AF chains along c axis are weakly 
coupled. However, for a > 0.3 the quantum corrections 
are already smaller than those in the G-AF phase, and 
further decrease with increasing a = J c / Jab, i.e. , when 
the FM interactions along the c axis get stronger. 

The stiffness constants D ab and D c shown in Fig. 
3(c) increase with increasing a in all cases, proving that 
the increasing exchange interactions cause magnons to 
harden. For weak inter actions along the c axis one 
finds D c /Dl D ~ 7^73 for G- and G-AF phase, and 
D ab /Dl D ~ v /q73 for A-AF phase, where D% D = 
2^/3JS is obtained in a G-AF phase with cubic sym- 
metry. We note that both D c and D ab increase when J c 



is increased in the G-AF phase. In contrast, the stiff- 
ness constants: D ab in the G-AF phase, and D c in the 
A-AF phase, are not influenced by the increasing FM 
interactions (J c and J ab , respectively), and one finds: 
D ab /Dl D = and D C /B% D = 1/V3 in the G-AF 

and A-AF phase. This illustrates a general rule - when 
the AF interactions increase, the spin stiffness increases 
both along the AF and FM bonds. In contrast, the spins 
on the FM bonds cannot fluctuate and thus the changes 
of FM interactions do not modify the spin stiffness along 
the AF bonds. 

The dependence of the quantum corrections and the 
stiffness constants on the interlayer coupling oc J c in the 
bilayer system is summarized in Fig. 4. As in the 3D 
case (Fig. 3), the increasing interlayer coupling gives a 
fast crossover from a system of two independent planes 
to a bilayer system with reduced values of S(S Z ) and 
SE/Eo [Fig. 4(a)]. It is quite remarkable that the ef- 
fect of finite and small J c is almost identical for the AF 
and FM interaction in the regime of \J C /J ab \ < 0.3, and 
the quantum corrections decrease simply due to the di- 
mensional crossover. Only at \J C /J a b\ > 0.3 the G-AF 
and G-AF phase start to differ: the quantum correc- 
tions pass through a minimum and increase again due 
to the AF interlayer coupling (J c > 0) which favors 
singlet states on these bonds in the G-AF phase, while 
they steadily decrease when the FM interlayer coupling 
{J c < 0) gets stronger in the G-AF phase. Similar as 
in the 3D case, D c /d\ ~ \faj2> in the regime of small 
a = \J C \/J ab < 0.3, where £>jj = 2\/3JS, found for the 



G-AF phase at J c = 2J ab 



is used as a unit. 

,11 



In the C- 



AF phase the value of D ab = a/2/3-Dq is independent 

of J c , while it approaches D ab = in the G-AF phase 
when |J C | — > 2J ab . Finally, as expected, the quantum 
corrections in the A-AF bilayer phase increase with in- 
creasing AF interlayer coupling J c > 0, and become even 
larger than those in the G-AF phase for sufficiently large 
J c . One finds equal reduction of the order parameter 
S(S Z ) ~ 0.095 in both phases at \J C /J a b\ ^ 2.46. 

The present study shows that the spin quantum correc- 
tions due to spin fluctuations are small both in the A-AF 
phase of LaMnOs, and in the G-AF phases of LaVOs 
and YVO3. Indeed, large magnetic moments <~ 3.87 
measured in LaMnOs (see Ref. 11) are almost perfectly 
reproduced by the present calculation which gives the 
moment (M) ~ 2(i B (S z ) ~ 3.91^ B for spin S = 2 and 
for the experimental ratio of \J a b\/J c = 1-43. In contrast, 
the strong reduction of the order parameter in LaVOa (to 
~ 1.3yU_B 12 ) and YVO3 (to <~ 1.0/x.b at finite temperature 
T ~ 85 K 13 ) cannot be explained by rather weak quan- 
tum spin effects, reducing the order parameter only down 
to - 1.89/is for spin S = 1 and | J c \/J ab ~ 1 at T = 0. 
This result indicates that the quantum effects due to or- 
bital fluctuations are stronger and dominate the behavior 
of the ti g systems with degenerate orbitals. 27 

Summarizing, we have shown that the quantum cor- 
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rcctions to the order parameter and to the ground state 
energy decrease systematically with the increasing num- 
ber and strength of FM bonds, from G-AF through C-AF 
to A-AF phase, both in the 3D cubic lattice and in the bi- 
layer system. This shows that the AF interactions dom- 
inate in the non-cubic antiferromagnets. Indeed, when 
the interactions are AF along only one or two cubic di- 
rections, their increase causes the increase of the spin 
stiffness and of the quantum fluctuations along all three 
cubic directions. 
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